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Although f{R) modified gravity models can be made to satisfy solar system and cosmological 
constraints, it has been shown that they have the serious drawback of the nonexistence of stars with 
strong gravitational fields. In this paper, we discuss whether or not higher curvature corrections 
can remedy the nonexistence consistently. The following problems are shown to arise as the costs 
one must pay for the f{R) models that allow for neutrons stars: (i) the leading correction must 
be fine-tuned to have the typical energy scale jj, < 10~^^ GeV, which essentially comes from the 
free fall time of a relativistic star; (ii) the leading correction must be further fine-tuned so that it 
is not given by the quadratic curvature term. The second problem is caused because there appears 
an intermediate curvature scale, and laboratory experiments of gravity will be under the influence 
of higher curvature corrections. Our analysis thus implies that it is a challenge to construct viable 
f{R) models without very careful and unnatural fine-tuning. 
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I. INTRODUCTION 



The origin of the current accelerated expansion of the 
Universe [1] is one of the biggest mystery in cosmology. 
The accelerated expansion may be driven by some un- 
known energy-momentum component. A more intrigu- 
ing possibility is that the acceleration could be due to 
long distance modification of gravity. A simple class of 
modified gravity theories can be constructed by general- 
izing the Einstein-Hilbert Lagrangian to some function of 
the Ricci scalar, f{R) 01. Various models of f{R) grav- 
ity have been proposedfl, 0] , but inappropriate choices 
of the function readily cause unwanted instability la 
or gross violation of solar system constraints [1, 0, 
The troubles arise due to an extra propagating scalar 
degree of freedom, and hence viable f{R) models must 
be constructed in such a way that the dynamics of this 
scalar field is carefully controlled. This is in principle 
possible, and indeed f{R) theories can be made to sat- 
isfy solar system and laboratory tests by invoking the 
chameleon mechanism 

i, Eli- The key ingredient 
of the chameleon mechanism is the density-dependent 
mass of the scalar field; it mediates a short-range force in 
high density environments such as the solar interior and 
vicinity. (The actual mechanism to hide the chameleon 
field is slightly more involved [Q] .) Concrete examples of 
"chameleon f(R)" are found in ^ El, El, El (see also 
Refs. [llE3,MEi,[13)- They are the only known exam- 
ples of viable f{R) models that exhibit no problems and 
no pathologies in the weak gravity regime [2ll [22. [23|.^ 
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^ The model of [l2j | , which belongs to a different class of the mod- 
els [T3I . [Tsl , is hardly distinguishable from ACDM cosmology 
because of the very strong experimental constraints [23] ■ 



However, the potentially viable models of [Tj, Ej, El 
turn out to have a serious drawback in the strong gravity 
regime. That is, a deep (but not diverging) gravitational 
potential drives the effective scalar degree of freedom to 
a curvature singularity. This problem was first pointed 
out by Appleby and Battye in a cosmological setting 1241 
and then discussed by Frolov in a general context [25|. 
In the previous paper [11, we have studied relativistic 
stars in f{R) gravity and shown explicitly that stars with 
strong gravitational fields develop curvature singularities 
and hence are prohibited. The critical value of the poten- 
tial is typically given by |$| ~ 0.1, implying problematic 
nonexistence of neutron stars in the models of [l3,E3,El- 

In this paper, we continue our program of study- 
ing strong gravity aspects of f{R) gravity, and discuss 
whether or not higher curvature corrections to the orig- 
inal models can resolve the singularity problem. This is 
done again by constructing relativistic star solutions. A 
higher curvature correction changes the structure of the 
effective potential for the scalar degree of freedom around 
the singularity [H, [27I. W e consider a modified version 
of Starobinsky's f{R) [ll, adding a correction term pro- 
portional to i?™ {m > 2). We also check whether or not 
the chameleon mechanism works to pass local gravita- 
tional tests in this modified f{R) model. Although we 
focus on the specific model, our result will hold in the 
other similar models of this class. 



This paper is organized as follows: In the next section, 
we describe the field equations of f{R) modified gravity 
in terms of a scalar-tensor theory. Then, in Sec. Illli we 
define the specific theory we consider. Our numerical 
results are presented in Sec. |TVl In Sec. El we discuss 
local tests of gravity in the f{R) model and point out 
the problem associated with the high energy correction 
term. We draw our conclusions in Sec. IVII 
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II. f(R) GRAVITY AS A SCALAR-TENSOR 
THEORY 



A. Field equations 

The action we consider has the form of 

fiR) 



S = d x^—g 



(1) 



where f{R) is a function of the Ricci scalar R, and is 
the Lagrangian of matter fields. Variation with respect 
to metric leads to the field equations^ 
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where fj^ := df/dR and T^^ := 
The trace of Eq. reduces to 



-2SCJSg>''^ + g^,C„ 



^{2f-fnR). 



(3) 



We now introduce an effective scalar degree of free- 
dom, which sometimes is dubbed "scalaron," by defining 
X ■= fn- Inverting this relation, the Ricci scalar can be 
expressed in terms of x: R — Q{x)- this way Eqs. ^ 
and ^ are equivalently rewritten as [28] 



XG' 



SttG, 



2r^dy 

3 dx' 



x'vix)S,'X^) 

(5) 



where the potential V is given by 



^(x) ^ ixQix) - fiQix))] , 



(6) 



and dV/dx = [2/(Q(x)) - xQix)] /(2x')- 

Eqs. (jl]) and ^ are equivalent to the Jordan frame 
equations of motion in the Brans-Dicke theory with lo — 
plus a potential V{x)- One can move to the Einstein 
frame by performing the conformal transformation g^^ — 
X9fj.i' with X = exp{^yl6^^G/3(j)), where is the canonical 
scalar field. The potential for (j> is then given by V{x{(l>))- 
However, we do not work in the Einstein frame in the 
following discussion. 



FIG. 1: The potential V. The inset shows the structure 
around the de-Sitter minimum. The potential of the original 
model (without i?™ term) is shown by a blue line for purpose 
of comparison. Parameters are given by X — 2, n = 1, m — 2, 
and £ = 5 X 10~^. The point X = 1 corresponds to a curvature 
singularity in the original model, but the i?™ term pushes the 
curvature singularity toward infinity, x ~ oo. 



radial profile of x through Eq. ([5]), it is useful to note 
that the equation can be written as 



fx 

dr^ 



2dx 
r dr 



dU ^ 
dx 



where 



dU 
dx 



2x^ dV 
~d^ 



(7) 



(8) 



and !F = {8ttG/3)T. Here we have ignored the effect of 
the metric for simplicity. (Later we will solve the full set 
of the field equations numerically.) Now, by identifying 
r as a time coordinate, Eq. ([7]) can be regarded as the 
equation of motion in classical mechanics. One can un- 
derstand the radial profile of x intuitively as the motion 
of a particle in the potential U under the time-dependent 
force !F (and the frictional force corresponding to the sec- 
ond term in the left hand side). The mechanical analogy 
is particularly useful to comprehend the essential point of 
the nonexistence statement for relativistic stars in /(i?) 
gravity [26l |. 



III. ADDING HIGHER CURVATURE 
CORRECTIONS TO f{R) GRAVITY 



B. Classical mechanical analogy 

We are going to investigate static, spherically symmet- 
ric stellar solutions in the above system. To study the 



^ In this paper, we focus on the metric approach rather than the 
Palatini one. 



In the previous paper [2(| we studied the strong grav- 
ity aspect of Starobinsky's f{R) theory described by 
/(i?) = i?-|-Ai?o[(l+^V^o)""-l] fll- There we showed 
that stars with strong gravitational fields (e.g., neutron 
stars) cannot exist in this model. We argued that this 
statement applies to the other similar models [3, [3 as 
well. This problem arises due to the dynamics of the ef- 
fective scalar degree of freedom, x, in the high curvature 
regime. Therefore, the problem may be cured by adding 
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FIG. 2; The effective potential U. The inset sliows tlie struc- 
ture around tlie de-Sitter extremum. Tlie effective potential of 
the original model (without i?™ term) is shown by a blue line 
for purpose of comparison. Parameters are given by A = 2, 
n = 1, m = 2, and e = 5 x iO~*. The dangerous curvature 
singularity is pushed toward x = oo by the i?*" term. 



higher curvature corrections that modify the structure of 
the potential around the large R region, as already noted 
in the original reference [l^ and later discussed in (27j . 

In general, higher curvature corrections may be writ- 
ten as a2R^ + asR^ + ■ ■ ■ , and so the most natural choice 
of the leading order term will be R'^/iJ?.^ The R^ term 
may be responsible for inflation in the early Universe if u 
is set to be an inflationary scale (e.g., ji ~ 10^^ GeV) [29j . 
but in this paper we do not restrict the mass scales of the 
curvature correction and assume that such parameters in 
the high energy correction terms take rather arbitrary 
values. If the coefficient of the R^ term is highly sup- 
pressed for some reason, then the leading correction will 
be the form of R^ /fi'^.^ To make the model simple but 
general enough, let us consider a function 



fiR) = R + XRo 



^1 -1 
Rl 



R" 



2(m-l) ' 



(9) 



where n{> 0),A(> 0), Ro{> 0), and m (> 2) are pa- 
rameters. The present Hubble scale is basically given 
by Hq ^ 0{Ro). We define a dimensionless param- 
eter e := Ro/lJ'^ and assume that e ^ 1 since the 
last term in Eq.([5]) is the high energy correction. At 
sufficiently low energies we have no cosmological con- 
stant, /(i?) ~ R, while for Rq <t: R <t: we find 
f{R) ~ R- XRo + AE^"+V^^" + ■ • ■ ■ At very high 
energies, i? » /x^, the last term dominates. 



^ Higher order corrections naturally include terms like R^^R^"', 
but in this paper we focus on the /(i?)-type modified gravity 
and hence simply assume that the corrections are also given by 
a function of the Ricci scalar. 

* Note that the higher curvature correction is given by the quartic 
terms in type II superstring theory. 



A de Sitter solution, R = Ri = constant, minimizes 
the potential V{x), and hence is found by solving the 
algebraic equation 



2fiRl)-R,fR\R^R,^0. 



(10) 



We may define the effective "cosmological constant" as 
Aeff := i?i/4. 

The scalar field x is written in terms of R as 



X 



1 - 2nX 



R 

Rn 



Rl 



(11) 



In the original model without the i?™ correction, a cur- 
vature singularity R = oo corresponds to a finite x 
(x = 1) and this is very close to the de Sitter mini- 
mum, xi — x(^i)- However, as is clear from Eq. 
the dangerous curvature singularity now corresponds to 
X = oo, and hence one may expect that this model is 
safe. A typical form of the potential V{x) is shown in 
Fig. [T] The eflFective potential U{x) is also plotted in 
Fig. [2l A straightforward calculation shows V oc 
for R^ Therefore, V — > const, as _R ^ oo for m = 2, 
while 1/ — > in the same limit for to > 3. Similarly, we 
have 



3— Pi-R+{m-2)- 
dx 



2(m-l) 



(i?»i?o). (12) 



From this we see that in the to = 2 case the effective 
potential U becomes steeper as the curvature increases, 
leading finally to dU/dx —oo as i? ^ -|-oo. For to > 3, 



U has a minimum at R 
R +00. 



fi and dU/dx — > +oo 



IV. RELATIVISTIC STARS IN f{R) GRAVITY 
WITH HIGH ENERGY CORRECTIONS 

We now investigate static and spherically symmetric 
stars with constant densities (i.e., a generalization of the 
Schwarzschild interior solution) in the model defined by 
Eq. (IHj). We shall work along the lines of the previous 
paper [2a]. The basic equations are found there and are 
replicated in Appendix Stars in f{R) gravity have 
been studied also in Ref. .30]. 

Given a density po and the central values of the pres- 
sure Pc and the scalar field Xc (or, equivalently, the cen- 
tral curvature Rc), we can integrate Eqs. (|A3|) - (|A6p nu- 
merically from the regular center r = to the surface 
of the star, r = TZ, defined by p{TV) = 0. (The bound- 
ary condition at the center is also given in Appendix [Xl) 
Then, imposing the continuity of the metric functions 
N{r) and B{r), the scalar field x, and its derivative dx/dr 
at the stellar surface, we integrate the vacuum field equa- 
tions (|A4[) - (jA6|) to find the exterior geometry. We are 
looking for a solution such that it is asymptotic to de 
Sitter with Aes — Ri/4: (and hence x ~* Xi)- For fixed 
Pq and Pc, we can find the desired solution (if it exists) 
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FIG. 3: Plots of the Ricci scalar R{r) for different e. Param- 
eters are given byA = 2,n = l,m = 2. The energy density is 
AivGpo ~ lO^Acft and the central pressure is pc = O.Spo- Solid 
(dashed) lines correspond to the region inside (outside) the 
star. These examples typically give GM/TZ ~ 0.25 - 0.26. 
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FIG. 4: Plots of xi^) for different e. Parameters are the same 
as those in Fig. |31 Solid (dashed) lines correspond to the 
region inside (outside) the star. 



FIG. 5: Plots of the metric functions for e = 5 x 10~^. Solid 
(dashed) lines correspond to the region inside (outside) the 
star. 
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FIG. 6: Plots of R{r) for e = 5 x 10"^°. Parameters are the 
same as those in Fig. [3] The upper (red) line is a plot for the 
solution with the central curvature Rc — 0.1185x87rGpo. This 
solution exhibits the behavior of falling toward the singularity. 
The lower (blue) line indicates a "overshooting" solution with 
Rc — 0.1184 X SnGpo. The inset emphasizes the overshooting 
behavior. 



by carefully tuning the initial value Xc = Xcrit- In the 
mechanical analogy, this solution corresponds to the sit- 
uation where the particle starts at rest and reaches the 
top of the potential (x = Xi) in the limit of r — > cx). The 
particle overshoots the top of the potential for Xc < Xcrit, 
while it turns around before it reaches the top and falls 
into the singularity for Xc > Xcrit- 

In the previous paper j26| we showed that Xcrit becomes 
larger as the gravitational potential of the star increases, 
getting eventually at Xs, above which the slope of the 
potential dU /dx is greater than the force term and so 
the particle cannot climb up the potential. Thus, if the 
gravitational potential is larger than a certain value, the 
desired solution described above does not exist, and we 
only have two types of singular solutions, i.e., a "falling- 
down" type and a "overshooting" one. 



Bearing the above mechanical picture in mind, let 
us move on to the case with the i?™ correction. We 
have carried out numerical calculations for various val- 
ues of e = i?o//i^ with fixing the other model parame- 
ters as A = 2, n = 1, and m — 2. The profiles of x 
and R for regular, asymptotically de Sitter solutions are 
shown in Figs. [3] and 21 and the metric functions for the 
e = 5 X 1G~^ case are plotted in Fig. [S] In these plots 
the energy density and the central pressure are given re- 
spectively by A-kGpq — IG^Acff and Pc ~ 0.3po, leading 
to the gravitational potential as large as GM/TZ ^ 0.25, 
where G := G/xc and M := Anp^Tl^/'d. (The gravita- 
tional potential is controlled by the ratio Pc/pa-) Stars 
with such large potentials are prohibited in the original 
model without the high energy correction. 
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For e < 5 X 10~ , however, we find that the regular 
solution ceases to exist. To see this, we show the behavior 
of the Ricci scalar for e = 5 x 10^^" in Fig. [6l Taking 
Rc = 0.1185 X SttGpo 1.185 x 87rG'(po - 3pc)], X goes 
toward x = oo (i? = oo). In this case, the pressure is an 
increasing function of r away from the center, and there 
is a maximal circumferential radius corresponding to an 
infinite proper distance from the center, as can be seen 
from Fig. [7l It is not an asymptotically flat space, but a 
"cylindrical" shape space with a singularity at the end. 
Taking a slightly smaller value, Rc = 0.1184 x SttGpq, x 
then overshoots the top of the potential and rolls down to 
the left. During this rolling-down phase the Kretschmann 
scalar, R^^^pcrR^^"'"^ , diverges. Only these two cases are 
realized, both of which are unphysical. This is the same 
situation one encounters in the model without the high 
energy correction [2^. In order for the higher curvature 
term to come to the rescue, it must be sufficiently large. 

The minimum value of e (or the maximum value of /i^) 
that allows for relativistic stars depends on the energy 
density. To explore the bound on /x^, we have performed 
numerical calculations for different values of po ranging 
from Po = 10'^Aeff/(47rG) to po = 10^Aeg/i4.TTG).^ As 
an example, the behavior of x for po = 10^Acff/(47rG') is 
shown in Fig. [S] From our numerical results it is con- 
firmed that the minimum value of e is inversely propor- 
tional to the energy density (Fig. [5]), and we roughly have 
the bound e > IO^^Rq/ (SirGpo) in order for stars with 
strong gravitational fields (~ 0.25) to exist. ^ This condi- 
tion gives 

p^ <ax SnGpo, a - 0(10^). (13) 

Taking Po - pnuci lO^-^g/cm^ - lO^^ Qgy'^ (nuclear 
density), one arrives at /i < 10^^^ GeV. This result itself 
is not surprising because SnGpnuci is a natural scale as- 
sociated with neutron stars. However, purely from a the- 
oretical point of view, this provides an unnaturally small 
energy scale. Obviously, the R^/p^ term with such small 
p cannot be relevant to inflation in the early Universe. 

We have also done numerical calculations to construct 
stellar solutions with GM/TZ ^ 0.26 in the m — 3 model, 
and obtained essentially the same result: for sufficiently 
large e we can find regular, asymptotically de Sitter so- 
lutions, while for e smaller than a certain value we only 
have two classes of unphysical solutions. Note here that 
for m > 3 the structure of the (effective) potential near 
the curvature singularity is quite different from that of 
the m = 2 model. Nevertheless, the solution correspond- 
ing to xC*") moving toward right shows an unphysical na- 
ture: the pressure is an increasing function of r away from 



^ To mimic a neutron star, the density must be SirGpo ~ lO^^A^ff . 
However, it is difficult to implement such an extremely high den- 
sity contrast in our numerical computations. 

® For given e there is a maximal gravitational potential. Therefore, 
the minimum value of e will be different depending on how large 
gravitational potentials one needs. 
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FIG. 7: If the force J-' is too weak to defeat the potential 
slope dU/dx inside a star and consequently x{r) and R{r) 
grow without turning back, the pressure shows an unphysical 
behavior. The coordinate choice is not good in this case, 
as can be seen most clearly from the behavior of the metric 
component Qrr- Plots are for \ — 2, n = 1, AnGpo = lO^Acff, 
Pc = 0.3po, and Rc = 0.12 x SnGpo- 



the central region and grr — > cxd at finite r. To allow for 
relativistic stars with GM/TZ ^ 0.26, it is required that 
/i2 < a' X SttGpo where a' - 0(10). 

Before closing this section, let us comment on the be- 
havior of the metric for asymptotically de Sitter stellar 
solutions. A numerical fitting leads to the approximate 
expression for the metric outside stars: 

A^ - A^oo (^1 - 2ci^ - yAeffr^^ , (14) 

B ~ i_2c3--^Aeffr2, (15) 
r 3 

where C2 — C4 ~ 1.0 irrespective of e (> 5 x 10~^), while 
ci and C3 are slightly different for different e. For e = 
5 X 10"^ one finds Ci ~ 0.29 and C3 ~ 0.24. This gives 
the post-Newtonian parameter 7 ~ C3/C1 ~ 0.81. For 
e = 5x 10"^ one has a ~ 0.30, C3 ~ 0.21, and 7 ~ 0.72. 
These results imply that the chameleon mechanism does 
not work in the above examples. This is because we are 
considering a vacuum exterior. Note that taking into 
account the effect of surrounding media does not remedy 
the nonexistence of relativistic stars for small e: to avoid 
falling down toward large R, x inevitably overshoots the 
top of the potential also in the presence of exterior matter 
as we have T < outside the star. 
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FIG. 8: Plots of x{r) for AttGpo = IO^'AcH and = O.Spo- 
The model parameters are given by A = 2, n = 1, and m = 2. 
In the £ = 5 X 10~^^ case, the desired solution is obtained for 
Rc = 0.2245 X SnGpo- However, in the case of e = 5 x 10~^^ 
only singular solutions are found. Two examples are shown: 
one is for = 0.1185 x SnGpo (falling rapidly down to i? = 
oo) and the other is for Rc — 0.1184 x 8nGpo (a overshooting 
solution) . 



M -9 " 

o 



-10 



-11 



-12 1 



5 6 7 8 9 

logjQ[47rGpo/Acff] 



FIG. 9; The minimum value of e as a function of po- Points 
indicate numerical results, showing the scaling relation. 



V. APPEARANCE OF AN INTERMEDIATE 
SCALE AND ANOTHER FINE-TUNING 

Assuming that the "UV scale" is given by fi^ ~ 
87rGpnuci, one may naively expect that higher curvature 
corrections have no impact on local tests of gravity since 
relevant densities are much smaller. However, as we show 
below, this expectation is not true. The purpose of this 
section is to point out a new problem brought by the 
higher curvature correction. 

To discuss the behavior of gravity in laboratories, it is 



convenient to define the effective potential for the x field 

as 

^ = I mQix)) - xQix)] + ^T, (16) 

where the energy-momentum tensor of matter is in- 
cluded. We consider the regime Rq R <^ fi^. In this 
regime Eq. (jlip can be written as 



X 



1 - 2nX 



Rq 



2n+l 



R 



(17) 



R ) 

At the minimum of the effective potential V^cff one finds 



dx 







R K, ^-kGp, 



(18) 



where p sa —T is the energy density of nonrelativistic 
matter. The mass of the excitation of the x field around 
the minimum is given by = (PVc(i/dx^\Rad,TvGp- 
Thus, the Compton wavelength A^,- 
puted as 



^ can be com- 



fcl 

R 







fR\ 




'r 





,(19) 



where ki :— 6n{2n + 1)A and ^2 '■— 3m(m — 1). 

To evaluate the Compton wavelength, it is important 
to note that there is a critical curvature scale defined by 



l/{2n+m) 



(20) 



and for i? <C i?* (respectively, R ^ R^,) the first (re- 
spectively, second) term in Eq. is much greater 
than the other. Eq. (pO)) gives an intermediate curva- 
ture scale Ro <^ R* fi^, which implies that the high 
energy correction term comes into play in determining 
the Compton wavelength of x o,i o. much lower scale than 
expected. In terms of energy densities, the intermediate 
scale R^,/{8nG) may be written as 



l/(2n+m) 



(21) 



where poE ~ 10 ^"g/cm^ and puv ^ Pnuci ~ 
lO^^g/cm"^. For instance, putting n = 1 and m = 2 
yields < lO^^^g/cm"^. Laboratory experiments are 
usually done at densities much higher than this! 

Although the original model is made to satisfy so- 
lar system and laboratory tests, the intermediate scale 
brought by a high energy correction term can destroy 
its success. Indeed, in the above example (n — 1 and 
m — 2), the mass of x is independent of local energy den- 
sities for p and so the chameleon mechanism does 
not work in laboratories. Since x h^-s a gravitational- 
strength coupling and the Compton wavelength is eval- 



uated as Av 



P' 



"1 > 10^ cm, the 



m 



2 model is ruled 



out by the fifth force constraint [31|, [32 
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As seen from Eqs. p9|) and (|20p. the intermediate scale 
and its consequences are sensitive to the exphcit form 
of the higher curvature correction. For example, in the 
n = 1 and m = 3 case one obtains < 10~^^ g/cm'^, and 
for densities higher than this the Compton wavelength is 
found to be Xy^ > 0.1 mm x (p/1 g • cm^"^)^/^. This typi- 
cally gives the marginal scale tested by laboratory exper- 
iments of gravity. Thus, determining whether or not a 
given high energy correction satisfies local tests requires 
a more careful study, which is beyond the scope of the 
present paper. We just emphasize here that "high energy 
corrections" play a crucial role above the intermediate 
curvature scale, and the B? I ^J? correction, which seems 
to appear in natural circumstances, is clearly inconsis- 
tent with laboratory tests if one chooses the parameter 
/i^ so that the theory evades the nonexistence statement 
of neutron stars. 

Finally, let us comment on cosmology with the i?™ 
term. In the matter-dominated era, we have an esti- 
mate m^/iJ^ ~ {fjL^ / H'^)™-~^ > 1, where it is assumed 
that the matter energy density is much greater than p,. 
This implies that the excitation of x is suppressed, ren- 
dering the field safe for cosmology. Before the time of 
matter-radiation equality the energy density of nonrel- 
ativistic matter is given by pm = rpr, where is the 
energy density of radiation and r := a/aeq <SC 1. Since 
R ^ SnGpin and ^ SnGpr/S, one ends up with 
m^/H^ ^ r^"'+^(/i^/i?^)™~^. Even at nucleosynthesis, 
the ratio m^/H^ is enhanced by the factor r~™+^ ex- 
cept for m — 2. Thus, we can approximately reproduce 
standard cosmology for < p^ in models with m > 3. 

VI. SUMMARY AND CONCLUSIONS 



energy correction completely destroys the success of the 
original f{R) model that passes local tests of gravity. 
If the quadratic correction is suppressed relative to the 
other higher curvature terms, possibly this is not always 
the case. However, it might be thought of as a prob- 
lem that gravity in the intermediate curvature regime 
(i?o ^ i? <C p^) is so sensitive to the explicit form of 
UV correction terms. In this sense we need another fine- 
tuning of the high energy correction. 

To conclude, although there is still a very small room 
for a possible construction of viable f{R) models that 
evade local tests of gravity and allow for stars with 
strong gravitational fields, very careful and unnatural 
fine-tuning is required for the model construction, leaving 
challenges for f{R) modified gravity. 
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APPENDIX A: SPHERICALLY SYMMETRIC 
STARS IN f{R) GRAVITY 

In this appendix we summarize the basic equations 
for constructing spherically symmetric stellar solutions 
in /(i?) gravity 



In this paper, we have tried to resolve the singularity 
problem arising in the strong gravity regime of otherwise 
viable f{R) theories. Adding a higher curvature correc- 
tion in the form of we have studied stars 
with strong gravitational fields which were prohibited in 
the original models. Our numerical results have shown 
that the scale p^ cannot be taken to be as large as an in- 
flationary energy scale nor a natural UV cutoff scale like 
(SttG)^^. Rather, p'^ < 0{87rGp), where p is the stellar 
density and hence is taken to be a nuclear density, is re- 
quired in order to remedy the nonexistence of relativistic 
stars. This provides a "high" energy scale as small as 
p < 10~^^ GeV. This is the flrst fine-tuning required for 
the high energy correction. 

In contrast to the naive expectation, the high energy 
corrections come into play at an intermediate curvature 
scale in determining the mass of x's excitation around 
the minimum of the effective potential. If the leading 
correction is given by the quadratic curvature term, the 

intermediate scale is i?* ~ (i?Q"'^^/i^)^^^^"'''^\ and the 
corresponding energy density is p■^, ^ lO^^^g/cm'^ for 
n = 1. For densities higher than this, the Compton 
wavelength of x is P~^ ~ 10^ cm. Therefore, the high 



1. Basic equations 

We take the ansatz of a spherically symmetric and 
static metric: 

els'" = -N{r)dt^ + ^ + ^2 (de^ + sin^ Od^^) . (Al) 
B[r) 

The energy-momentum tensor of matter fields is given by 

T^"^ = diag(-p,p,p,p). (A2) 

From the energy-momentum conservation, V^T^" — 0, 
we obtain 

p' + ^ip + p)^0- (A3) 

Here and hereafter a prime denotes differentiation with 
respect to r. The (tt) and (rr) components of the field 



8 



equations (|4|) yield, respectively, 



in the power series of r as 



{-1 + B + rB') 



-B 



'l^^)x' . (A4) 



The equation of motion for x [Eq- ©] gives 
B 



X 



N' B' 
8nG 



X 



3 (-p + 3p) + — — . (A6) 

We do not integrate the angular components of the field 
equations. Instead, we use them to check the accuracy 
of our numerical results, because those are derived from 
other equations via the Bianchi identity. 

If the energy density is constant inside the star, p — pQ, 
Eq. (jA3|) immediately gives 



7V(r) 



Po+Pc 



Po+Pir) 



(A7) 



In the main text we only consider constant density stars 
for simplicity. 



2. Boundary conditions 

Let us study the boundary conditions at the center of 
a star. Assuming the regularity, we expand the variables 



N{r) = 1 + N2r^ + B{r) = 1 + Bar^ + 



X{r) = Xc ( 1 + + 



P{r)=Pc + Y'^ P[^)=Pc + Y^ + 



(A8) 



where Xc, Pc and Pc are the central values of the scalar 
field, the energy density and the pressure, respectively. 
Note that using the scaling freedom of the t coordinate, 
we set A^(0) = 1. From Eqs. (|X4l) - (|M)) . we obtain 

3^2 = -87rGpc-XcVc-5C2, (A9) 
B2 + 2N2 - SirGp, - XcVc - 2C2, (AlO) 



8ttG, 



2x1 



3^2 - — (-Pc + 3p,) + ^T/;,,, (All) 



where G := G/xc, K := V{xc), and F^c = dV/dx\x=xa- 
These three equations are rearranged to give 



B2 - -^(2p, + 3p,)-f K-%^K,.,(A12) 

8nG Y 
N2 = _(2pe + 3p,)-^K-^^^xc, (A13) 
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G2 = ^i-Pc + 3pc) + ^Vx^ 
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2x| 
9 



(A14) 



Then, p2 is derived from the conservation equation: 

P2+N2{pc+Pc)^0- (A15) 

The Ricci scalar is given hy R = Rc + 0{r^) with Rc = 
-6(^2 + N2) near r = 0. 
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